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THEORETICAL ANALYSIS OF THE LONGITUDINAL BEHAVIOR OF AN 
AUTOMATICALLY CONTROLLED SUPERSONIC INTERCEPTOR 
DURING THE ATTACK PHASE 1 

By Ordway B. Gates, Jr., and C. H. Woodlixg 


SUMMARY 

A theoretical analysis ha* been made of the longi- 
tudinal behavior of an automatically controlled 
super son ic interceptor during the attack phase. 
The control system used to control the interceptor’s 
flight path was one in which a pitching velocity was 
com manded in proportion to the longitudinal tracking 
error. Throughout the investigation the assumption 
is made that the target is flying on a straight-line 
path. 

Factors considered in this investigation included 
e fleets oj control -system parameters , effects of limita- 
tions on control deflection and rate of control de- 
flection , effects oj initial tracking errors, effects of 
nonlinear variations in drag and lift with angle of 
attack and Mach number, effects of nonlinear varia- 
tions in pitching moment with angle of attack , 
effect of variations in interceptor forward velocity, 
and the effect of a normal acceleration limiter on 
the system performance. 

The control system considered in this investigation 
was found to give acceptable control oj the inter- 
ceptor's flight path during attack runs against a 
no n ma newer ing target. 

The inclusion oj a nonlinear variation of drag 
and lift with angle oj attack and Mach member 
resulted in relatively large variations in the inter- 
ceptor forward velocity during the attack phase. 
However, the effects of velocity changes on the 
overall responses during the attack phase were 
considerably reduced when a signal proportional 

' Su^K’i’snlos ivmitly declassified NAC'A Research Memorandum L54K0N 


to the change in forward velocity was fed back to 
the elevator servo. 

INTRODUCTION 

In ail interceptor research program engaged in 
at the Langley Aeronautical Laboratory of the 
National Advisory Committee for Aeronautics, 
one of the purposes was to evaluate t lie tracking 
performance of a supersonic interceptor equipped 
with various types of automatic, control systems. 
The present paper is concerned with an analysis 
of that phase of the problem wherein t ho inter- 
ceptor's radar locks on, with an initial vertical 
tracking error, to a bomber (lying at a constant 
velocity; only maneuvers of the interceptor in 
the vertical plane are required to carry out the 
interception. The results obtained from analysis 
of this longitudinal phase of the general tracking 
problem are intended to provide information 
which will be useful in the synthesis of a satis- 
factory longitudinal control system for the inter- 
ceptor being studied. The interceptor considered 
in this investigation is similar to that analyzed in 
reference 1, which has a notched delta wing of 
aspect ratio 3.2 and 55° sweepback of the leading 
edge. 

For this investigation the interceptor is assumed 
to be Hying initially in level flight at a Mach 
number of 2.2 at an altitude of 50,000 feet, and 
the target is flying in level flight toward the 
interceptor at a Mach number of 1.4, at various 
altitudes above 50,000 feet. No consideration 

IV Ordway B. (tab's, Jr., and C, II. Woodlfng, 1955. 
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was given to the effects of altitude changes on 
the interception problem discussed in this paper. 

The guidance equations presented in tins paper 
are for a lead-collision type of navigation. 

The results of this investigation are presented, 
for the most part, in the form of interceptor and 
kinematic responses subsequent to radar lock-on 
and were computed on the Reeves Electronic 
Analog ( 'omputer (REAC). 


SYMBOLS 


. . Drag 

(\> trim drag coefficient, — ^ 


f/> per radian 

1 a da 

(\ trim lift coefficient, 
( =~- L y per radian 

L “ da 1 

pornuiian 


Lift 

qS 


(' m pitching-moment coefficient, 

Pitching moment 

qSc 

< = d< 67 ’ per r,ulial1 

d i>r 

< ,= ^ > per radian 

' Hr) 

r m per radian 

da 

> per radian 

a _ aC 
a 21- 

r -*,“ d ^’ pcrr,,di “" 

c mean aerodynamic chord, ft 

]) differential operator, — 

q acceleration clue to gravity, ft/sec 2 

J Y moment of inertia about 1 stability axis, 

slug-ft 2 

K t racking-loop gain constant, radians/sec/ 

radian 

K r rate of pitch-feedback gain, radians/sec*/ 

radians/sec 

K s elevator-servo gain constant, radians/radian 

/see 


M 

M 


d/.VL.S 


m 

n 

9 

It 


S 

t 

to 


u 

u' 

V 

a 



6 

a 


T 


T J 


n 


Mach number 

predicted miss distance, measured positive 
from interceptor to target, ft 
component of M along the instantaneous 
line of sight, positive when target is ahead 
of rockets at predicted time of impact, ft 
component of M perpendicular to the in- 
stantaneous line of sight, positive when 
target is below rockets at predicted time 
of impact, ft 
mass of airplane, slugs 
normal acceleration, g units 
dynamic pressure, lb/sq ft 
distance* from interceptor to target along 
line* of sight, measured positive from 
interceptor to target, ft 
wing ami, sq ft 
time, se*e* 

time* of (light of interceptor from instan- 
taneous position to firing point, sec 
change* in forward velocity, ft/soe 

. u 

relative change m forward. velocity, y. 


forward velocity, ft /sec 

angle of attack, radians unless otherwise 
specified 

flight-path angle a), radians unless 

otherwise specified 

elevator deflection, radians unless otherwise 
specified 


-K 

output of filter, radians 

error in interceptor’s (light path at any 

M r 

given instant , e, - , g “ 

angle* of pitch, radians unless otherwise* 
specifieel 

angle* between interceptor A body axis and 
radar line of sight, positive when line of 
sight is above body axis, radians unless 
otherwise specified 

time of flight of interceptor’s rockets from 
firing point to predicted point of contact 


with target, sec 
filter time* constant, see* 

[‘levator servosystem time constant, see 
angular velocity of lim* e)f sight, (12= <j +0), 
radians/see; positive* when line of sight is 
rotating upward 
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Subscripts: 

I interceptor 

/ input 

L limit 

0 initial value 

R rocket 

xx steady state 

T target 

ANALYSIS 


with the target to the predicted point of impact. 
The geometry of the attack problem is shown dia- 
grammatieally in figure 1. Generally, the vector 
equation which must be satisfied is: 


R+ V T (t a+r) - Vit G + ( Vj+ V*) r+M ( 1 ) 

The components of this vector equation along and 
normal to the instantaneous line of sight are: 


DERIVATION OF GUIDANCE EQUATIONS 

The type of navigation or interception con- 
sidered in this investigation is lead collision; that 
is, the interceptor endeavors to fly a constant 
flight path such that at only one point on the 
path the rockets of the interceptor may be fired 
and a hit obtained on the target. The rockets, 
subsequent to firing, fly a constant bearing course 


R-\~ T' V(f g~\~ T ) CDs (c r-\-0 — 7 r) ^ 

— [ V it <7+ ( V I h 1 r) t] COS (<t-\-ol) 4- M ls 

y 

1 r(t g~\-t ) sin (cr-f-0 y T ) 

= 1 1 jt <v+(l / + \ r) t | sin (<j +a) + M XL s y 


( 2 ) 


The target flight-path angle y T is taken as zero 
when the target is in level flight going away from 
the interceptor and is taken as w when in level 



Fkicke 1. Geometry of lead collision navigation used in present investigation for Mu g = 0. 
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flight coming toward the interceptor. The eqna- 
tions may be rewritten in terms of the range, rate 
of change of range, and angular velocity of the 
line of sight as: 

1\ + It ( to + t ) — 1 rt cos (<t+ a) = M Ls 
— RQ{t (j T t) — I rt sin (a T<*) = A/y ls 

where 


R=Yr cos (<r + 0 — 7 y-) — F/ cos (<T~ra) 



r Sill 

(cr T a) — I 7 - 

sin (< 

rT 0 — 

7 t ) 

quant it it 

*s It, 

u, 0, and a 

are d( 

:'finod 

as: 

7?= 

/fo + 





cr — 


| a(H=tr,r f 

[(0- 

-e)dt 


0 = 

0() + 

| Mt 




a — 







In practice, R, cr, and a would be available from 
the radar; however, in the Analog solution of the 
problem these quantities were obtained from 
actual integrations. 

Certain simplifying assumptions were made in 
this investigation. The angles (<r+0) and (cr + a) 
were assumed to be small enough so that the 
cosines and sines of these angles are equal to unity 
and to the angle in radians, respectively. For 
these assumptions 

R ~ 1 t [cos 77 t T(^ - 1 - 0) sin 7 y] — T / 
li$t~Yi (o’Ta)— Fr [(ff-f-0) <‘Os y T — sin y T ] 
and the guidance equations are 

R + R (ta+T)~^ V ~M ls ^ 

r ( 4 ) 

— r t (<r+a) J 

The solution of equations (4) is accomplished by 
computing continuously from the first of equa- 
tions (4) the value of (fi ; + r) necessary to make 
A/.vl.s — 0 and then, from the second of equations 



(4), computing for this value of (ton-r) the value 
of M XLl : which will exist at (fp+r) seconds sub- 
se(pumt to the instantaneous time. The time at 
which t a is computed to be zero is the firing point 
for the interceptor’s rockets. The time of flight 
of the interceptor's rockets is r, and throughout 
this investigation is assumed to be 1.5 seconds. 
The command to the control system is based, on 
the error € 7 which exists at any time in the inter- 
ceptor's flight path, which for this investigation 
is approximated by the expression 

M.wls /~\ 

Cr ~ Vjt 0 +{Vi+V R )r W 


The validity of the foregoing assumptions in the 
guidance equations were checked by a digital 
solution on the Bell Telephone Laboratories 
X-66744 relay computer at Langley where the 
exact guidance equations were used. This com- 
parison is discussed in a subsequent section. 


DISCUSSION OF FLIGHT-PATH CONTROL SYSTEM 


The block diagram of the overall system is 
presented in figure 2. Briefly, tin* computed 
(plant it v e 7 is filtered, amplified, and list'd as the 
comma id to a pitch-rate command system. The 
dynami *s of the filter and elevator servo are 
represented by simple first-order lag networks of 


1 


the form — and 

1 -j- T fiJ 

transfei function 


1 


v respectively. The 


1 4“ r 


T S W 

is assumed to he repre- 


sentatb e of a low-pass filter wliicli in practice 
would 1 e necessary to attenuate the high-frequency 
radar l oise present, in the computed command 
signal 7 . However, no attempt was made to 
include noise in the present, investigation. The 
dynamics of the interceptor wore obtained from 
the linearized equations of longitudinal motion. 
For eeitain cases, these equations were modified 
to include specific nonlinearities. All equations 
used ii the analysis and the analog schematic 
diagram are presented in the appendix. The 
interceptor parameters and other constants used 
in the analysis are presented in table I. The 
interceptor stability derivatives and mass param- 
eters were obtained from unpublished data and 
the results of reference 1 . 
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Target 



Figure 2. Block diagram of the longitudinal tracking system used in present investigation. 


TABLE I RESULTS AND DISCUSSION 

SELECTION OF SYSTEM GAIN CONSTANTS AND EFFECTS OF 
VARIOUS AERODYNAMIC PARAMETERS ON 
ATTACK PERFORMANCE 

System gain constants. — An investigation was 
first made to determine values of the gain con- 
stants K and K r for which the at tack performance 
of the interceptor would be reasonably sat isfactory 
as a starting point for the general study. The 
servo gain constant K, is taken as unity through- 
out this analysis. Examination of the block 
diagram (fig. 2) indicates that the forward -loop 
gain is KK S and the feedback gain is K r K s ; there- 
fore, the assumption that K s equals unity imposes 
no restrictions on the system gain constants or 
performance. From a theoretical analysis of the 
open-loop frequency response yjt y and several 
preliminary runs on the REAt 1 , acceptable values 
of K and K r were found to he 3.0 and 0.375, 
respectively. 

For these values of K and K r tracking runs were 
computed for 7? 0 = 60,000 feet, and <7 (J — 7.5° 

and 15° and are presented in figure 3. For these 
runs the interceptor motions were computed from 
the linearized equations of longitudinal motions 
presented in the appendix. The responses shown 
in figure 3 include the predicted miss distance 
normal to t lie line of sight M NLS) the interceptor 


STABILITY DERIVATIVES AND MASS 
CHARACTERISTICS OF INTERCEPTOR AND 
OTHER CONSTANTS USED IN INVESTIGATION 

Altitude, ft 50, 000 

Density p, slugs/eu ft . .. 0. 0003022 

r /n , ft /sec _ 2, 140 


M _ 

m, slugs 

ly, slug-ft 2 

q , lb/sq ft 

7, ft 

,S, sci ft 

C„, q: per radian _ 

C tn u , per radian- 
cy, per radian _ 

( , , j)cr radian. _ 
C Lat per radian _ _ 
C i ) a , per radian __ 

C„ 

C L 

C m , per radian__ 

6 r 

Cl j per radian __ 

U/e, ft/sec* 

I t, ft /sec 

t s , see 

ry, see 

t, see 

/C, ft 

On, radians 

radians 


2 . 2 

770. 4 

2. 08 X 1 0 5 
820 
If) 
401 

_ - -2. 84 
_ . - 0. 50 
- -0.28 

0 . 00 

_ - 2. 29 

0. (50 

0. 027 

_ 0. 070 
- 0. 295 

. 0. 105 

_ 2 , 000 
.. 1 , 300 
.. . 0. 03 

0. 00 
1.5 

. 00, 000 
.. 0. 033 
0.033 
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Tracking responses were computed for A”— 3.0, 
K r — 0.375, fl„=(K),()0() feet, and * 0 =7 .5° and 15° 
for tlie cases where equations (6) and (7) were 
used for C D and C La . The results are presented in 
figures 4(a) and 4(1)). Also shown in these figures 
is the case for which AC d is assumed to vary 
lineaily with Aa, that is, AC d — 0.156Aa. For 
<7 o= 7.5° (fig. 4(a)), the change in interceptor for- 
ward velocity is seen to be 0.071 A ,o for the linear 
case, ().ldV/,o when A( f 1} varies nonlinearly with 
Aa , and 0.14 V r /t() when A( varies nonlinearly with 
Ac* and AM and C L varies wit h Ma<*h number. For 
the linear cast* M XLS — — 50 feet, and for t he non- 
linear eases 4 / jV ls ~ —50 feet . As pointed out previ- 
ously, the nonzero values of M#ijs are due to the 
change in forward velocity. r Flie stunt' general 
trends are noted for <r 0 =15° (fig. 4(b)), but the 
changes in forward velocity for this value of <r 0 
are much greater for each drag condition investi- 
gated than were encountered for <7 {) = 7.5°. Also, 
the values of A /. v ls at the end of the run are 
larger for this value of <j () than for <r n =7.5°. 

Effect of feedback proportional to change in 
forward velocity. — A possible means of eliminating, 


Linear A Cq, linear 

Nonlinear A Cq [a Cq - /(Aa)j , linear C ^ 
Nonlinear A Cq [A Cq : /(Aa, AAU], 



J 1 I I X 1 i 



/, sec 

(a) Co — 7.5°. 

Fun re 4. -Effort of nonlinear drug and lift on interceptor 
5:t74r»;i <io 2 


or at least reducing, the value of M NLS at the as- 
sumed firing time is to feed back a signal to the 
elevator servo proportional to the change in for- 
ward velocity, such that a positive pitching 
moment is produced for a decrease in forward 
velocity. Time histories are presented in figures 
5(a) and 5(h) for <70 = 7.5° and 15°. For these 4 
runs, C D varies nonlinearly with Aa- and AM and 
a feedback gain of 0.12 is utilized. For each 
value 4 of <r 0 the predict eel value of Mxi^ is seen 
to be appreciably reduced. This type of feedback 
requires a bias error in the (light path in order for 
the 4 interceptor to fly a constant flight path; but, 
on the basis of the predicted value of M NLS pre- 
sentod in figure 5, this bias appears to be small. 
For comparison, runs m which the change in for- 
ward velocity was assumed to he zero are also 
shown in these figures. Despite the fact that the 
velocity changed considerably during these runs, 
there appears to he no appreciable difference 
between the cases which included the feedback 
proportional to u' in which forward velocity was 
allowed to vary and in which velocity changes 





/, sec 

(l>) ffo 15°. 

attack performance. K -3.0; K r 0.375; /A = 00,000 feet. 
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(a) Assumed variation of (■ m with or. 

(b) Time history. 

Figure 0 .— Effect of .nonlinear pitching-moment coefficient 
with angle of attack on interceptor attack performance. 


increase the static sensitivity between pitching 
velocity (or normal acceleration) arid elevator 
deflection. This effect is reflected in the larger 
normal acceleration for the nonlinear case. 

GENERAL EFFECT OF CONTROL-SYSTEM PARAMETERS AND 
INITIAL CONDITIONS ON ATTACK PERFORMANCE 

Tracking runs were computed and the results 
indicate the general effect on the system per- 
formance of variations in the gain constants and 
other system parameters. The comparisons are 
presented in the form of time histories and in the 
form of summary plots which indicate the effect 
of the various quantities on factors such as rise 
time, response time, and maximum overshoot with 
respect to the miss distance M X ls, and maximum 
normal acceleration encountered during the run. 


In this discussion, the rise time is defined as the 
time required for M NLS to reach initially the point 
of zero miss and the response time is defined as the 
time for M XLS to reach and remain less than 30 
feet. 

Effect of the tracking loop gain K. — The effect 
of the gain const ant K on the system performance 
is shown in figure 7 for values of <r 0 equal lo 7.5° 
and 15°. The results are presented for K— 3, 5, 
and 9. For these (rases the initial range was 
60,000 feet and the feedback gain K r was equal to 
0.375. The effect of K is seen to be essentially 
the same for both values of < 7 0 . As the gain is 
increased the system becomes more and more 
oscillatory and for K— 9.0 the M XL s response is 
probably unsatisfactory for both values of cr 0 . 
The effects of K with respect to the rise time, the 
response time and overshoot of the M\ls transient, 
and the maximum normal acceleration encountered 
are summarized in figure 7(c) for a () — 7.5° and 
15°. As K is increased, there is only a slight 
variation in rise time and peak acceleration, which 
is due, to a large extent, to tin* fact that the 
(devator reaches its maximum deflection of —20° 
very quickly and remains at that deflection for a 
t ime dependent upon K and <r 0 . As K is increased, 
the response time hauls first to decrease and then 
to become large as K is further increased. The 
response time for K 9 and cr„ 15° is not shown, 
since for this combination of K and <x 0 , M sls 
never reaches the condition whore it remains less 
than 30 feet. For both values of co, the over- 
shoot in A /.van increases progressively with in- 
creases in K. 

Effect of rate feedback gain K r . For K— 3.0, 
< 7 0 — 7.5°, and 7£ n = 60,000 feet, results are pre- 
sented in figure 8(a) for values of K r — 0, 0.20, 
0.375, 0.60, and 1.0. For K r — 0 the responses are 
seen to be rather oscillatory. As K r is increased 
the system becomes more stable, hut the over- 
shoot. in Mnlh i s seen 1° increase with increases 
in I\ r . These effects are due to the fact that, as 
the interceptor is stabilized, its response to control 
inputs becomes slower and results in an increase 
in the M XLS overshoot. A summary of the effects 
of K r on the A/ V /,n transient is presented in figure 
8(b). Inclusion of the pitch-rate feedback tends 
to reduce the response time and overshoot of 
A v ls and, in addition, to eliminate the oscillatory 
condition which exists for K r — 0; but for values of 
K t greater than 0.375, the overshoot is larger than 
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for K t 0. There is seen to I >e a relatively negli- 
gible effect of K r on the rise time because, as 
mentioned previously, 5 V is at its limit of —20° at 
tilt' beginning of the run and the rate feedback 
is ineffective during the early portion of the run. 

Effect of initial error in ov — For i£=3.0, K T — 
0.375, and Ii 0 ~ 00,000 feet, several runs were made 
to evaluate the ability of the system to score a 
hit as the initial error in a is increased. The 
value of a actually reflects the altitude difference 
and the horizontal distance between the inter- 
ceptor and target ; at the beginning of the run the 
alt it mb' difference is given by the expression 
I{ u sin(<r o + 0o). Results are presented in figure 9 

for <r„ = 2°, 7.5°, 10°, 15°, and 20°. As <r 0 is 
increased, the time which elapses between reaching 
the correct flight path (Afxis — ft) and the firing 
point becomes less and less, until for a {] — 20° the 
interceptor is unable to reduce d/ v/ ,. s to zero. 
The maximum allowable initial value of cr 0 is a 
function of the range and the maximum accelera- 
tion that the airplane can pull. In the present 


system for 5, f/ , = — 20°, is approximately %. 

Presented in figure 10 is a plot, of the initial 
range I\ } against the maximum initial <r 0 for which 
the present, system could score a hit. This 
curve v as obtained from a simultaneous solution 
of the equations 

VrJ ( OS [y(t)\dt + (Vi + VR)T cos [t (to)] 

— V T COS y T (ta+ t) = lit) cos O o +0<>) 

y (*) 

V i | ^iu \y{t)\dt 4- ( l / V R ) r sin |q(f a) 1 

J n 

— Vt si i 7 = sin (<r o + 0n) 

for the case where y T — tt. Equations (8) relate 
the horizontal and vertical distances traveled by 
the interceptor, rockets, and target to the hori- 
zontal and vertical distances which exist between 
the inlerceptor and target at If these 

equations were satisfied, a hit would be obtained. 
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(e) Summary plots. 
Ficit : rk 7. — Concluded. 


The function y{t) was calculated from the longi- 
tudinal equations of motion for a step input on 
8 e = — 20°, which is taken as the maximum value 
of 8 C throughout these calculations. Also pre- 
sented in figure 10 is the variation of /i 0 with 


oo which was calculated on the assumption that 
y(t) — 7(0) + 7/. The value of 7 used in this 
expression is the steady-state 7 due to an elevator 
deflection of —20°, and it. was assumed that the 
airframe attained this output immediately upon 
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(a) Tiim* history; 7.5°. 

Kum kk S. Ktfeet of pitch rate feedback on interceptor 
attack performance. K 3.0; 1>0,00() feet. 


application of control. The results obtained for 
I his simplified approach to the problem are seen 
to he substantially the same as those obtained 
when the more exact approach is used. The 
curve indicates that for AV- : 60,000 feet the largest 
a 0 which can be used is approximately 22°. A s 
a cheek, runs were made for a large value of K 
(in order to keep <5 e — (5 f )z, for the entire run) 
and the results indicated that <r 0 — 20° was about 
the largest value which could he tolerated. 

In order to obtain a hit for A\ ( = 60,000 feet and 
00 — 22 °, the required time of flight of the inter- 
ceptor from its initial position to the assumed 
firing point is seen from figure 10 to he 16.8 
seconds. This explains the difference between 
the maximum allowable <r 0 indicated by figure 
10 and the value obtained from the REAC 
results. For the KEAC runs, the simplifications 
made in the guidance equations eliminated the 
dependence of (f G +r) on <r 0 ; and for the initial 
range A’,,— 60,000 feet, t a is equal to 14.8 seconds. 
Actually, t G for the RE AC runs was slightly less 
than 14.8 seconds, since li {) was closer to 59,000 
feet, due to voltage limitations on the RE AC. 





(b) Summary plots. 

Figure 8. — Concluded. 

A check was also made for /f 0 = 30,000 feet and 
the ret ult also agreed well with the curve of 
figure 1 1 0. 

Effect of filter time constant r f . — The effect of 
T f is illustrated in figures 11(a) and 11(b) for 
T f — 0, 0.30, 0.60, and 1.20 seconds. For these 
runs /v — 3.0, K r = 0.375, <r 0 =2° and 7.5°, and 
/?o=60,000 feet. As r f is increased from 0 to 
1.20 seconds, the initial response in M NLS is seen 
to become progressively slower, which is due to 
the increased lag between the initial command 
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2 

7.5 

10 

15 

20 





Fkjijrk 9. -Effect of initial values of radar elevation 
angle a on interceptor attack performance. K 3.U; 
K r -=0.375; = < >0,000 feet. 



Eictjre 10. — Combinations of /fo, <ru, and to for which a 
hit may be obtained in this investigation. (a 0 — 0 o— 1 .ft° ; 
7 T = 7T-) 


e y and the response oi the elevator motion; there 
is also seen to be an appreciable overshoot in 
the M N ls response for r f — 1 .20 seconds for the 
same reason. The effects of r f on the At sls 
are summarized in figure 11(c) for a 0 ~2° and 
7.5°. 

The values of a {] chosen to illustrate the effect 
of r f have no special significance, but were used 
only on the basis that the results obtained for 


these values of <j {) were typical of the results 
obtained for all values of a (i up to the maximum 
allowable value for AV- - 00,000 feet. The same 
statement may he made concerning Ihe values 
of (To used in the subsequent sections. 

Effect of servo time constant r s . The effect of 
the servo time constant r s is demonstrated in 
figure 12 for r,-(U)3, 0.10, 0.20, and 0.80 second. 
The eases presented are for /C=3.0, A^ r = 0.37o, 
a o ” 2 ° , and /A, —00,000 feet. As r v is increased 
from 0.03 second to 0.30 second, the most im- 
portant effect is that the overshoot in the Af SLS 
response is seen to become larger, but at the as- 
sumed firing point, the miss distance is zero in 
either case. The effects of r v on the AI XIS response 
are summarized in figure 12(b). 

Effect of control-surface limitations. —The 
effects of limitations on the rate of control-surface 
deflection and the magnitude of the surface 
deflection were investigated briefly for K— 3.0, 
K r = 0.375, and AV= <>0,000 feet. All of the 
results presented in this paper up to this point 
were obtained for the condition where was 
limited to ±120°/sec and 8 e to ±20°. The effect 
of reducing the maximum value of 5, to ±()0°/see 
and ±30°/sec for this tracking run may bo seen 
from figure 13(a). For these runs <r {) - 10°. The 
limiting control deflection was kept at ±20°. 
The effect' on this run of reducing (5,)/, was to 
reduce slightly the peak normal acceleration as 
(6 e ) L is reduced from ± 120°/see lo ±30°/sec. 
This reduction in acceleration causes a slight 
increase in tin* rise and response time of Afy LSy 
but, it may be concluded that, at least for this 
ease, the effects of limitations on the rate of 
control deflection were small. However, for 
cases where higher values of K would be required 
(for example, a maneuvering target) the effects 
of rate limitation would probably be much more 
important and should ho investigated thoroughly. 
The effects of limiting 8 e on the A f XLS responses 
are summarized in figure 13(b). 

Several runs were made for the ease of (d r ) L 
— ±120°/sec and (Sj L reduced from ±20° to 
±10°; the results are presented in figure 14. 
For these cases <r 0 is equal to 5°. The general 
effect of reducing the control deflection limits 
from 20° to 10° is to reduce the peak normal 
acceleration of the interceptor, and hence the 
rise and response time of tin* Af XLS response. 
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r/, sec 
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CD 

01 

-o 


CO 


t , sec 

(a) Time* history; a 0 2°. 

Kioi kk II. KftVi't of fil tor-tii no constant r/ on inlem*ptor attack performance. 


t , sec 

(1)) Time history; ao 7.5°. 

K 3.0; K r 0.375; /iV 00,000 foot. 



The maximum a u for which M.\ls <*an ho reduced 
to zero is directly dependent on the limit im- 
posed on 5 r . 

Effect on tracking of limiting the command#,, 
Since the result obtained from filtering and 
amplifying the error signal e y is used as the com- 
mand to a pitch-rate command system, this 
command should be limited if it is desired to 
limit the interceptor’s normal acceleration from 
aerodynamic, pilot comfort, structural, or other 
considerations. The normal acceleration re- 
sponse to 6i is given by 



lienee the steady-state normal acceleration can 
be limited to any desired value by limiting the 
input command B t . For the interceptor being 
considered 

( 6 \ ... L 

W„ W9)..+K r K, 4.85+A, 


an<l, in general, 



Therof »n>, 

(fff)i=(w)t(4.8S+A r )-^ (!<)) 

If Bi limited by this expression, it should he 
pointe 1 out that only the steady-state u is being 
limited, by use of this expression. The effective- 
ness of limiting the output transient acceleration 
by this method depends primarily on the response 
characteristics of (///#,)• If the system gains arc* 
cliosei: to give a fast response of n to 0* with little 
or no overshoot, this means of limiting n should 
be satisfactory. The results presented in figure 15 
afford a comparison of the cases for which there is 
no //-limiter and for the case where <j is limited 
by equation (10). For these cases K— 3.0, K r 
=-0.375, 7^ = (>0,000 feet, <r (J = 15° and (n) L ~rn/. 
The lirframe steady-state* normal-acceleration 
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r s , sec 

0.03 

.10 

.20 

.30 




(a) Time history; cr„ — 2°. 

Kiouhk 1 2. -Effect of servo time constant r 3 on inter- 
ceptor attack performance, /v -3.0; K r = 0.375; — 

00.000 feet. 



response to the limiting value of 5 t is approxi- 
mately 4.0//; lienee, the unlimited // ease and the 
ease of (n) L — 5// should oscillate about tin* same 
value when 5 t —(5 ( .) L . As can be seen, the peak n 
for the unlimited ease is roughly N//. When (0 f )n 
is limited by (‘([nation (10), the peak // response is 
reduced roughly to 0.5//. The rise time of 4 f VLS 
is seen to increase slightly when is limited, but 
the response time is less than that for the unlimited 
ease. 

Digital check on validity of simplified guidance 
equations. In order to cheek the validity of the 
small-angle assumptions made in the guidance 
(‘([nations (see previous section entitled “Deriva- 
tion of Guidance Equations”), a solution, using 
the exact guidance of kinematic* equations, was 
obtained from the Bell Computer for comparison 
with the REAC tracking solutions which utilized 
the simplified equations. This comparison is 
presented in figure 10 for the case of o- 0 = 7.5°, 
TO, 000 foot , /v = 3.(), and 1C = 0.375, and the 
agreement is seen to he excellent. 



(b) Summary plots. 
Figure 1 2. — Concluded. 
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, deg /sec 



/, sec 


(a) Time history: <r ( >=lO°. 

Fiot kk 13. — Kffoct of limiting rate of elevator deflection 

d r oil interceptor attack performance. JC=3.0; K r = 

0. 375. tf (J = 60, 000 feet; (8 e )/,= ± 20. 

CONCLUSIONS 

The following conclusions were reached from a 
theoretical investigation of the longitudinal track- 
ing behavior of an automatically controlled inter- 
ceptor against a nonmaneuvering target : 

1. The control system considered in this inves- 
tigation (i.e., command on rate of pitch propor- 
tional to longitudinal tracking error) was found 
to give acceptable control of the interceptor's 
flight path during attack runs against a non- 
maneuvering target. 

2. The inclusion of a nonlinear variation of drag 
and lift with angle of attack and Mach number 
resulted in relatively large variations in the inter- 
ceptor forward velocity during the attack runs. 

3. The changes in forward velocity computed 
for the runs in this investigation, although rather 
large, had a relatively small effect on the overall 
responses when a signal proportional to the change 
in forward velocity was fed hack to tin* (‘levator 
servo. 

4. ( Consideration of a nonlinear variation of 
pitching moment with angle of attack which 
tended toward static instability at high angles of 
attack indicated that its primary effect was to 





60 120 
(8*) , deg/sec 

(») Summary plots. 

Kkjurk 13. — Concluded. 


increase the magnitudes of the interceptor’s 
motions during the tracking runs. 

5. The general effect of increasing the tracking 
gain K was to destabilize the tracking loop. 

G. Increases in the rate-feedback gain K r tended 
to stabilize the interceptor's longitudinal short- 
period oscillation, hut had a destabilizing (‘fleet on 
the tracking loop. 

7. The maximum initial angularity between the 
interceptor's flight path and radar line of sight for 
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(8»)i, deg 



t, sec 


Fi'M kk 1 4. KfFoet of limiting elevator deflection (5 ,)l 
on interceptor attack performance. K— 3.0; K r -~- 0.375; 
AV=6 0,000 feet ; <7,1 = 5°; («„)/,= ± !20°/sec. 



0 4 8 12 

t, sec 

Fioi kk 15 — Comparison of unlimited 0, and 9 t limited to 
value miuired for «* 8 = 5. 



/, sec 


Fiut-RK 1(>. — -Comparison of interceptor attack perfor- 
mance for exact and simplified guidance kinematics. 
K = 3. ); A\ = 0.375; ff 0 = 60,000 feet; <x 0 = 7.5°. 

which a hit cun he obtained ran be well approxi- 
mated from the initial range and normal accelera- 
tion capabilities of the interceptor. 

8. increases in either the filter time constant 
r f or the servo time constant t s had an adverse 1 
effect on the attack performance because of the 
increased lag between the input command and 
the elevator motions. 

9. T ic effects of limitations on the rate of con- 
trol defection did not appear to be large in this 
investigation for the limiting rates considered. 

10. 'die general effect of limitations on the 
magnitude of the control deflection was to slow 
down the interceptor’s responses, and lienee the 
maxim im initial tracking error which can be 
tolerated decreases as the limits on elevator de- 
flect ioi are reduced. 

11. Limiting of the interceptor normal accelera- 
tion wi s achieved by limiting the input command 
on pi t < hing velocity, which effectively limits the 
intern* otor’s steady-state normal acceleration. 

Lano le y Research ('enter. 

National Aeronautics and Space Administration, 
La sOLKY Field, Va., October 22, U154- 



APPENDIX 

EQUATIONS USED IN ANALOG OF TRACKING PROBLEM AND ANALOG SCHEMATIC DIAGRAMS 


All equations used in this investigation are pre- 
sented in this appendix. The equations are pre- 
sented in both general form and with numerical 
substitutions for the parameters in table I. It 
will be noted that certain of the equations pre- 
sented in numerical form have been multiplied bv 
constants. These constants were used in order to 
adjust the REAC voltages to satisfactory levels in 
the analog procedure. 

Linearized airframe equations (7 />0 = 0): 


■■ 7 qSc • . s , 7 qSc . 

m 9 1 ~7 7 & 5 ^ '«« 9 f r t a 
v & V / o Jr - 1 /,o 1 y 

qSc , , qSc , ,,, , qSc , 

m a j a ~T ' rn . . t ~ i ■“ U i ' m§ i ~ O e 
ly * Y e Jy 


a = e-(\ 2 (\ <! f. u' 

“ rn\ ‘ml-,,. 


/ 1 / *7^ x / 


Y (Al) 


or - </S v' r <lS ft 

mV / o m\ /.o 


( 


. «« 






m V,,r ' * 


It will be noted that ( O', and 6/ appear 

in those equations rather than (\„ 6 , C Lh , and 8 e . 
The relations exist that 


8/ = -8 e 


O/ — — O,, 


This convention is adopted in order that the true 
physical phase relationships of the pitch-rate com- 
mand loop be obtained for positive gain constants. 


Control equations: 

. e 7 € f 

€/ T f T l 
B t—Ktf 

\t, r., / r. 

Kinematic equations (simplifiecl, 7 =tt) 

r 

1{= — V T — V T 

Ril ~ V / (cr -f- a) 1 t(& 

Yf VLS = [ r) 1 i/T (a -|-a) 

. Mnls 

V ilta • t) • t r t 

V, = V l . a + V,.„u' 


t, 


R /f„+J/A// 

a = a 0 -|- -^adi 

e- h// 


-J» 

■f 

<r = <r„+J(S>-e)r// 


/>’S2 

R 


(A2) 


(A3) 


For the parameters presented in table I, equa- 
tions (Al), (A2), and (A3) take the following form: 
Airframe equations: 

6 = - 0 . 1 84 50 - 0 . 0 1 8 1 9a - 1 0 . 382a -j~ 5 .4(195, ' 
a=0—O.45(L5a— 0.0303//' -f 0.03295/ 
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— 2()«'— 0.21 ■>»/.' +o.:«):«+o.r»79a— <).3<)3a 


('out rol ('([Hat ions: 

-t/~ 3.333*-,- -3.333*, 

6 j — 

O.OoS,/— 1 .(>(>70, — 1 X\iuK r 6 — 1 .(>(>76/ 

Kinematic equations: 

3000- R 

I (?+T= ; 

n 

12/?=, -16,320- 12 V, 

20/t > <2 — 42,800 ((7+a) + 27,200 (a+6) 
— 4-U vi6 .=4/i > lI (t 0 +T)+ 12,000 (( r+a) 

_~2M xls 

~ fy " — r/Oo-fr)— 3000 

-201 -42, 800-42, 800«' 

-A’ -00,000— J ihlt 

— 2a — 0.000— 2 iarit 


—0= -0.033- 



a — a 



-6) dt 


-1012 = 


-10 Ml 

R 


The a and 6 appearing in equations (A3) are 
total a and 0, whereas the a and 6 in equations 
(Al) are perturbations away from the trimmed 
condition. The solution was slowed down to the 
extent that 2 seconds of machine time was 
equivalent to 1 second of problem time. The 
scale factors used were: 100 volts =100 seconds, 
100 volts =1 radian, and 100 volts =00,000 feet. 

The analog schematics of equations (Al), (A2), 
and (A3) are presented in figures 17(a), 17(b), 
and 17(c), respectively. 
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SYMBOLS 



Multiplying servo 



integroting 

amplifier 


Potentiometer 


(a) 


(a) Airframe equations ( VI). 

Fuii KK 17. Analog schematic diagram of equations (Al), (A2). and (A3) and table of schematic symbols. 
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(b) Control equations (A2). (See fig. 17(a) for key.) 
Figure 17. — Continued. 
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